Semiclassical Description of Tunneling in Mixed Systems: 
the Case of the Annular Billiard 
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We study quantum-mechanical tunneling between symmetry-related pairs of regular phase space 
regions that are separated by a chaotic layer. We consider the annular billiard, and use scattering 
theory to relate the splitting of quasi-degenerate states quantized on the two regular regions to specific 
paths connecting them. The tunneling amplitudes involved are given a semiclassical interpretation 
by extending the billiard boundaries to complex space and generalizing specular reflection to complex 
rays. We give analytical expressions for the splittings, and show that the dominant contributions 
come from chaos- assisted paths that tunnel into and out of the chaotic layer. 
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Recently there has been a surge of interest in quantum 
systems whose classical counterparts exhibit a mixture 
of regular and chaotic motion. Such systems are of great 
importance since they comprise the majority of dynam- 
ical systems found in nature. One interesting feature is 
quantum-mechanical tunneling between regular regions 
in phase space that are separated by a chaotic layer. Nu- 
merical studies of mixed systems have shown that the 
splitting of quasi-doublets associated with such pairs of 
regular regions was much higher than could be explained 
by direct tunneling processes ■ This was attributed 

in to a suggested mechanism of chaos- assisted tun- 
neling, i.e. tunneling from one regular region into the 
chaotic sea, propagating "classically" to the other side, 
and tunneling out into the second regular region. Since 
a large part of the phase space is thus traversed via clas- 
sically allowed transitions, these paths were expected to 
considerably enhance the splitting. 

In this work we study tunneling from a scattering the- 
ory point of view, taking as a specific example the annular 
billiard proposed in || . This allows to write the splittings 
as a sum over paths in angular momentum space, which 
can include both classically allowed and tunneling tran- 
sitions. The tunneling amplitudes are evaluated semi- 
classically by continuing the billiard boundaries to com- 
plex configuration space, and generalizing the mechanism 
of specular reflection to the case of complex rays. We 
obtain analytical expressions for the splittings in terms 
of these transition amplitudes, and show that the domi- 
nant contributions arise from paths that tunnel into (and 
out of) the chaotic layer via intermediate angular mo- 
menta which lie on the boundary between the regular 
and chaotic regions. This approach yields values for the 
splitting which are in good agreement with exact results. 

The annular billiard consists of the space between two 
non-concentric circles of radius R and a < R, centered 
at (x, y) coordinates O = (0, 0) and O' = (—6, 0), respec- 
tively. Classically, a particle moves freely between spec- 
ular reflections on the bounding circles. We parametrize 
trajectories by their impact parameter L with respect to 



O and their direction of propagation 7. Trajectories of 
\L\ > a + S do not hit the inner circle, but rotate for- 
ever at constant L, while the phase space for \L\ < a + 5 
consists of a mixture of regular islands and chaotic lay- 
ers. Phase space plots can be found in Ref. |§. In this 
work we use parameters for which a single chaotic layer 
extends from L = a + 5 to L = — (a + 5). 

We apply the scattering approach to quantization first 
proposed in Q. Let us write the wave function in terms 
of incoming and outgoing cylindrical waves, 
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where H^ 1,2 ^(a;) denotes the Hankel functions of first and 
second kind, and k is the wave number. We consider 
the billiard as two back-to-back scattering systems; the 
first, "inner" system consists of incoming waves being 
reflected to outgoing waves by the exterior of the inner 
circle, whereas the "outer" system consists of outgoing 
waves being scattered to incoming waves by the interior 
of the outer circle. The two system are characterized 
by scattering matrices S^ l, °\k), respectively, that re- 
late the coefficient vectors a , (3 by (3 — S^'(k) a and 
a = S 1 (k) /3 . Requiring the two relations to be con- 
sistent results in the quantization condition 

det(£(fr) - 1) - , S(k) = S (I) (k)S (0) (k) . (1) 

Thus the billiard supports an eigenvalue whenever one of 
the eigenphases of S(k) equals an integer multiple of 2tt. 

We drop the wave number argument in the se- 
quel, is clearly just the diagonal matrix Snm = 
-E.£\kR)/Ki 2) (kR) S nm , while can be obtained 
by writing it as S n m = — H^ 2 ^ (ka)/ (ka) S nm in the 
primed coordinates, and then performing a coordinate 
change. Using the addition theorem for Bessel functions 
(see e.g. ||) this gives 

r (2) (/ca) 
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The structure of S can be seen in Fig. GL which shows 
|5 nm | for k = 100, a = 0.4 and S = 0.2. The inset gives 
an overview of one quadrant. One can see that S is al- 
most diagonal at angular momenta n,m > k(a + S). For 
n, m below this value most of the amplitude lies within 
the region of classically allowed transitions and is de- 
limited by caustics due to classical rainbow scattering. 
Closer inspection reveals that these ridges extend into 
the classically forbidden region as well, although they are 
exponentially suppressed there. The main figure shows 
a single row n = 70. Since n > k(a + S), the diagonal 
element is almost unimodular, while the non-diagonal el- 
ements are exponentially small. One finds a maximum 
at m w 63 that corresponds to the ridge seen in the in- 
set, while away from this maximum the matrix elements 
decay faster than exponentially. 

For our purpose it will prove sufficient to approximate 
the magnitudes \S nm \ — |<Snm|- We will now sketch 
this semiclassical derivation; a more detailed account 
will be given elsewhere 0. Starting from (^), we write 
Jm-e(kS) = e ™( m - e ) J e _ m (kS) and apply the Poisson 
summation formula to get 
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and consider only the /i = term (see below) . We replace 
the Bessel functions by their Sommerfeld representation 
J v (z) = l/(27r) J c d/y exp[iz cos j + iv (7 — 7r/2)], where 
the contour of integration C runs from — 7r+ioo to 7r+ioo. 
The contribution of the Hankel functions to non-diagonal 
elements of is asymptotically given by their De- 
bye approximation to be i exp[i 2fca(sin Q — Q cos Q)] , 
where Q = arccos(-^ / ka) . Snm can then be written in the 
form of an integral over an exponential whose argument 
is proportional to k. In the semiclassical limit k 3> 1 
this integration can be carried out using the saddle point 
approximation, giving 
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where the summation is over saddles p in the complex 
Q,7i,7/ space. These saddles are determined by 



, - kS (sin7/ )P - sin 7^) , 
-kS sin 7i p + ka cos Q p . 



(5b) 
(5c) 
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Solving (g) yields a number of saddle points in the com- 
plex plane. For most values of n, m the dominant con- 
tribution comes from a single saddle, whereas others are 
either negligible or cannot be reached by deformation of 
the contour of integration. This dominant saddle coa- 
lesces with its symmetry-related counterpart at values 
of 71, m which lie close to the ridges seen in the inset of 
Fig. 0. (The actual ridges correspond to points where 
7/ iP becomes real and equal to it/ 2 or 3tt/2). In princi- 
ple, near this caustic saddle point integration should be 
replaced by a uniform approximation, and also the Han- 
kel functions should be treated more carefully. However, 
a detailed treatment of this region is not required for 
the purpose of this paper. The full line in Fig. HJ shows 
the contribution of the dominant saddle to |5„m|, and 
one can see that there is excellent agreement over a wide 
range of angular momenta. 

For completeness, let us consider the fi ^ terms in 
@. Each integral in the summation can be written as a 
residue sum over the poles of B-P{ka)/ U^ 1} (ka). These 
residue contributions are interpreted in standard diffrac- 
tion theory as surface waves excited on the inner disc 
by rays of grazing incidence S. As can be seen from 
Fig. [l], at our parameter values they give rise to signif- 
icant corrections only over a small range of m. As the 
corresponding transitions will turn out to be of no inter- 
est in the present context, we defer a discussion of the 
residue terms to [||. 

The approximation given in Eqs. (|^), (^) can be inter- 
preted in terms of classical trajectories in complex con- 
figuration space, by extending the dynamics inside the 
billiard to complex coordinates. The extension of free 
flight to complex coordinates is trivial. In order to de- 
scribe the interaction with the inner circle, we define a 
complex circle of radius a around the point (xo,yo) as 
the surface for which x — xq = a cos (3, y — yo = asin/3, 
where (3 is allowed to be complex. Our motivation is that 
a wave function satisfying Dirichlet or Neumann bound- 
ary conditions on such a circle for real fi will also satisfy 
them for complex j3. Next, we note that through (al- 
most) every point on a circle one can pass two distinct 
trajectories with a given impact parameter relative to its 
center. Since reflection from a circle centered at the ori- 
gin preserves the impact parameter, we define specular 
reflection as the mapping from one of these trajectories 
to the other. 

It is easy to verify that any given initial and final im- 
pact parameters are linked by at least one complex trajec- 
tory. The conditions specifying the corresponding initial 
and final angles are then equivalent to Eqs. (||). Thus 7^ 
and 7^ iP have the meaning of the initial and final angles 
for which an incoming trajectory with impact parameter 
Li = n/k is reflected off the complex inner circle to the 
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impact parameter Lf = m/k. Moreover, ( paj ) is the gen- 
eralized specular reflection condition, while (|^,^) give 
the classical deflection function. Note that the impact 
parameter with respect to the center of the inner circle 
is given by V = acosQ. By (|5c[), we see that V can 
become complex even if Li and Lf are real. Further- 
more, unlike in J^], both the initial and final angles are 
generically complex. 

Finally, it is straightforward to show that the phase <& 
in (^) is given by the reduced action (not the length) of 
the ray, $ = - f* f dt [r(t)p r (t) + ip(t)L(t)}, where (r, ip), 
(j>r,L) are canonically conjugate polar coordinates and 
momenta, and that \1Z P \ is the corresponding reciprocal 
stability \d 2 $/dLidL f \/2<7r. We therefore find that (|) 
coincides with the sum over classical trajectories which 
constitutes the usual semiclassical approximation to scat- 
tering matrix elements (see Jul). 

We will now extract the level splitting. We first exam- 
ine the splittings in eigenphases, and will regard energy 
splittings later. Let | ± ) represent the doublet of eigen- 
vectors peaked at angular momenta ±n, 
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where the are exponentially small. We denote the 
corresponding eigenphases by 8^ and their splitting by 
59 n = \9+- 9-\. Using exp(iiV0±) = ( ± \S N \ ± ) we get 
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which describes the tunneling oscillations between n and 
— n, with a correction term < fc 2 max|K^| 2 . For 

we can therefore write 
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Eq. 



relates the splitting to a sum over all paths 



{Xi}lL 1 in matrix element space that go from Xi = —n 
to A^v = n. The simplest such paths involve only the 
direct tunneling amplitude S- n ^ n , which from Fig. ^ is 
exceedingly small. A much larger contribution comes 
from paths which tunnel only over small distances in an- 
gular momentum space and traverse the remaining dis- 
tance via classically allowed transitions. These are the 
chaos-assisted tunneling paths proposed in M. In order 
to deal with such contributions, we first need to consider 
the structure of the chaotic part of phase space. 

The quantum evolution of systems with purely chaotic 
classical counterparts is usually very well described by 
random matrix theory pT| . However, this is only appro- 
priate when phase space can be assumed to be completely 
structureless. In the case of a classically mixed system, 
this is an over-simplification: as the Lyapunov exponent 



vanishes smoothly at the interface between chaotic and 
regular phase space regions, there is always an interme- 
diate layer in which chaotic, but relatively stable motion 
gives rise to classical staying times much longer than the 
mean level density. Moreover, angular momentum re- 
mains a preferred basis well into the chaotic part of phase 
space. This corresponds to the observation made in [p||l2|] 
that classical regular structure can be quantum mechan- 
ically continued into the chaotic sea. Consequently, it 
is only the internal part 1 1 | < ^coe of the chaotic sea 
that is appropriately modelled by a random matrix en- 
semble, which in the present case is a circular orthogonal 
ensemble (COE) of size ~ 2Zcoe- 

It turns out that there are two types of chaos-assisted 
paths that contribute dominantly to the splitting: (I) 
paths (n, 7, —n) that tunnel from n directly into the 
chaotic region, propagate in some eigenstate I7) of the 
internal block, and finally tunnel to —n, and (II) paths 
(n, i, 7, — I, — n) that propagate from n to I7) and then 
to —n in two jumps via intermediate edge angular mo- 
menta I and — I, respectively. We will first consider paths 
of type (I). Let us denote the eigenphase of | 7) by 7 , 
and let 8 n = —i\ogS ntn . Summing over the possible 
dwell times at n and 7 gives contributions of the type 
Sn,^ S-y,- n / sm[(9 7 — 9 n )/2] to the splitting. Assuming 
that the S n ,-y are distributed independently of the # 7 , we 
can average over the COE and calculate the median of 
the type (I) splitting contributions |13| 
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where i* 2 = ( |5„ j7 | 2 ), which we approximate by 



< « E l g C =-i COE I S n , g I 2 - Similarly, the type (II) contri- 
butions can be estimated by 
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Note that S9 n will fluctuate strongly about the median 
values due to avoided crossings |^|,[|,[l4| . The relative im- 
portance of the contributions S9n' IJ ^ is system-specific; 



for the present system we find 69 { n H) /69 ( n> - 25. Fur- 
thermore, paths that include additional transitions out- 
side the chaotic block give negligible contributions. 

It remains to connect eigenphase and eigenvalue split- 
tings. From (0J|) we see that 59 n (k) is approximately 
constant over exponentially small ranges of k. It is then 
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where 9„ 
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simple to show that 5k n w | 39^ /dk 
is evaluated at 6 — 0. However, we expect Eqs. (f?],|]) to 
over-estimate the actual splittings: our model neglects 
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residual transport barriers inside the chaotic block, and 
the damping effect of imaginary parts of eigenphases is 
not accounted for. As these errors are expected to be 
independent of n, we can correct by an overall factor c 
that we extract from the numerical data. 

In Fig. we present a numerical test of our results. We 
plot median values of S6 n that were obtained by varying 
the outer radius R over 30 values between 1 and 1.3. 
Note that changing R leaves transition probabilities con- 
stant. The full circles represent the exact S6 n , as func- 
tion of n. Splittings smaller than ~ 10~ 14 (denoted by 
empty circles) could not be calculated directly due to 
the finite precision of the diagonalization procedure, but 
were obtained by calculating S N for large N and apply- 
ing Eq. (ph. Finally, the dashed line shows the prediction 
of Eqs. |f|) for Z COE = 50 and c = 0.1. While the n- 
dependence due to tunneling is reproduced very well, the 
coefficient due to transport across the chaotic sea could 
only be calculated up to an order of magnitude. The inset 
shows Sk n /S9 n evaluated from the quantization condition 
(0) around k = 100 for R = 1, compared to | 86^/dk I" 1 . 
We see that the correspondence is very good. 
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FIG. 1. Tunneling amplitudes |Snm| for n — 70, a = 0.4, 
S = 0.2 and k — 100 as a function of m, calculated exactly 
(dots) and semiclassically (full line). The inset shows one 
quadrant of \S\, with larger values corresponding to larger 
dots, in arbitrary units. Off-diagonal ridges at n, m > k(a + S) 
are exponentially enhanced in the plot. The dashed lines in- 
dicate k(a + S). 

FIG. 2. Median eigenphase splittings S8 n obtained by di- 
agonalization of S (full circles), using Eq. (|^) with TV = 2 13 
(empty circles), and as estimated from Eqs. (f?],[|) (dashed 
line) with c = 0.1 (see text). Inset: Sk n /S9 n (diamonds) 
and | d9rP/dk | _1 (dashed line) evaluated around k = 100. 
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